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We introduce an easily computable topological measure which locates the effective crossover be- 
tween segregation and integration in a modular network. Segregation corresponds to the degree 
of network modularity, while integration is expressed in terms of the algebraic connectivity of an 
associated hyper-graph. The rigorous treatment of the simplified case of cliques of equal size that 
are gradually rewired until they become completely merged, allows us to show that this topological 
crossover can be made to coincide with a dynamical crossover from cluster to global synchroniza- 
tion of a system of coupled phase oscillators. The dynamical crossover is signaled by a peak in 
the product of the measures of intra-cluster and global synchronization, which we propose as a 
dynamical measure of complexity. This quantity is much easier to compute than the entropy (of 
the average frequencies of the oscillators) , and displays a behavior which closely mimics that of the 
dynamical complexity index based on the latter. The proposed toplogical measure simultaneously 
provides information on the dynamical behavior, sheds light on the interplay between modularity 
vs total integration and shows how this affects the capability of the network to perform both local 
and distributed dynamical tasks. 

PACS: 89.75.Fb, 05.45.Xt, 89.70.Eg 



Many physical and biological systems (such as elec- 
tronic devices, communications networks, and the hu- 
man brain) face similar constraints as they interact with 
complex environments, and organize their structure and 
function along similar principles of resource allocation 
On the one hand, the need for fast and reliable re- 
sponses to changes in the environment naturally favors 
the emergence of segregated modules of specialized com- 
putation (e.g. sensory systems in the brain). On the 
other hand, interactions among modules become essen- 
tial when an information processing whose complexity 
exceeds the capacity of the single modules is required. 
For instance, perceptual systems in the brain need to 
bind information from different brain areas to produce a 
single coherent percept 0. Therefore, segregation into 
specialized modules and integration into global coherent 
activity present an inherent trade-off, and an appropri- 
ate balance between these two tendencies has been shown 
to be necessary for efficient functioning, particularly in 
neural systems Q. In fact, an exceedingly segregated 
or integrated functioning of the brain has been associ- 
ated with variouspathological conditions, e.g. autism or 
schizophrenia [3-[6 | , and epilepsy Q respectively. 

One straightforward way to study such a balance in 
complex systems is to represent them as dynamical net- 
works, endowing them with well-studied topological and 
dynamical properties (see [H, H[ for a review). For in- 
stance, Zhao et al. characterized systems of coupled 
phase oscillators in terms of a complexity index based on 
the entropy of the distribution of pairwise synchroniza- 
tion. Heterogeneous and modular networks were shown 
to be characterized by high complexity, for intermediate 



levels of modularity, in a regime marked by the forma- 
tion of dynamical clusters and the coordination between 
them. 

In this Letter, we provide an easily computable topo- 
logical measure for quantifying the balance between seg- 
regation and integration in a network. We propose that 
segregation can be understood in terms of a community 
structure (i.e. clusters of vertices densely connected to 
each other while less connected to vertices outside the 
community [Hjj), while integration can conveniently be 
expressed in terms of algebraic connectivity [ijj of the 
hyper-graph associated to the network. After introduc- 
ing our measure for a generic modular graph, we focus 
on the simplified case in which the network communi- 
ties are of equal size, and have the same number of inner 
and outer connections, and show analytically that there 
is a structure that maximizes the product of segregation 
and coordination measures. We then demonstrate that 
the dynamics emerging from such a specific configuration 
is associated with the coincidence of the two thresholds 
for cluster and complete synchronization in a network of 
interacting phase oscillators. 

We start by considering a generic undirected, un- 
weighted graph Q composed of N nodes and L links, 
partitioned into C communities, and characterized by an 
associated N x N adjacency matrix A. Moreover, for 
each community r (r = 1,...,C), we consider the num- 
ber of links connecting pairs of members of that com- 
munity (i^ n ), and the number of inter-community links 
(^oui)i i- e - the number of links connecting a member 
of that community with a member of another commu- 
nity. Based on the mentioned notation, we have L = 
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(E^=i £r m ) + I (E?=i Crt) ■ Usin S thc previous defini- 
tions, the standard modularity measure defined in [l3| . 

can be written as Q = Y? r =Mn/L ~ [Win + Ct ) / 2£] 2 } ■ 
Let us now define the hyper-graph Q* associated to 
Q, as the weighted directed C-clique in which each node 
corresponds to a community of Q, and the connection 
incident to node r from node s is weighted by -jr 1 , being 
i r s ut the number of links of Q that connects members of 
the community r with members of thc community s, and 
£ r in the number of inner links in the source community. 
The corresponding C x C Laplacian matrix C* = {C* s } is 
asymmetric, but can be written as the product C* = BC, 
where C is a symmetric zero row-sum matrix with off- 
diagonal elements C rs = —£ r s ut and diagonal ones C rr = 
tout = E.^r P'ouu and B = diag{l/C ■ • • , 1/ig}: 
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The spectrum of C* is real with non-negative values 
and because C* is zero row-sum, the smallest eigenvalue 
A J is zero, while > 0. The measure that we propose 
for the balance between integration and segregation is 
defined as follows: 



£ = QA5. 

Indeed, while Q is an inherent evaluation of the seg- 
regation factor of a graph, > quantifies the con- 
nectiveness of the hyper-graph, and therefore measures 
the extent to which different communities are bounded 
and interact. It should be noticed that both Q and 
are properly normalized in such a way that, even if the 
network links were associated to cohesive forces, the two 
quantities would be a-dimensional. The maximum of £ 
corresponds to a topology in which integration and seg- 
regation have the same weight. 

Let us then consider the case in which the C com- 
munities are cliques having equal size N c = N/C, and 
the number of intra-community links, P~ inl as well as 
the number of inter-community links, i r out , are the same 
for all communities. We then have i r in = £i n /C, and 
Pout = tout/C, such that L = C(£l n + \t r out ) = t m + \l out . 
Under these assumptions, the modularity can be reduced 
to the following expression, 



Q = C 
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On the other hand, we have that B = -f^-X and, from 
matrix identity, 




FIG. 1. (Color online) (a)-(c) Graph representations of (a) 
a network consisting of two cliques of size Nc = 16, and 
(b-c) the resulting networks after rewiring, (b) j max — QQ 
(£out/2L = 0.5 and Q = 0), and (c) j = 30 (£ ovt /2L = 0.25 
and Q = 0.25). Intra-community (inter-community) links 
are light gray (black colored). Blue (red) nodes are phase 
oscillators whose natural frequency is randomly distributed 
around 0.8 (0.2) within a range ±0.1. (d)-(f) Time evolution 
of the instantaneous frequencies of the iV = 32 oscillators af- 
ter solving Eq. ((4} with the networks of the left panels and 
din = d* n = 0.14. (d) d out = 0, (e) d out = d* out = 1.10, and 
(f) dout = 1.40. Dashed line in (d) marks the mean frequency 
between the two clusters, (g)-(h) R (squares) and (circles) 
(see text for definition) for: (g) the two clique network of (a) 
as a function of di„, (h) the 0- modularity network of (b) as a 
function of d out and for di n = 0.14. The vertical dashed line 
in (g) marks the critical coupling for cluster synchronization 
(d* n = 0.14), while in (h) it corresponds to d ou t — 1.4, well 
above the onset of global synchronization (d out = 1.10). Panel 
(i) reports the S and I indexes (see text for definitions) vs. 
the ratio £ ou t/(2L), after solution of Eq. Q with di n — 0.14 
and d ou t = 1.40. Oblique dashed line is the corresponding 
modularity index Q. 
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Let us now consider a particular protocol by means of 
which l out is varied from to 2L(C — 1)/C, the value 
at which the modularity Q is zero. We start from a 
fully segregated configuration (in which £ ou t = and 
£i n = L), and operate successive rewiring processes, in 
each of which an intra-community link from each commu- 
nity is deleted, and C inter-community links are formed 
by connecting those pairs of nodes (each one in differ- 
ent communities) having lost their intra-link. In this 
way, at the j-th rewiring, we have £i n = L — Cj, and 
P-out — C2j. The maximum number of rewiring steps 
until modularity fades out is, therefore, N C (N C — 1)(C — 
1)/(2C). It follows from combining Eqs. ([T]) and © 
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FIG. 2. (Color online) (a) Complexity index _E (circles) and S- 
I index (triangles) (see text for definitions) corresponding to 
the data reported in Fig. [Hi), (b) S ■ I curves as a function of 
the mixing ratio £ out /(2L) for C = 2 and for increasing values 
of d ut (in the legend), (c) Maxima of the S-I curves in (b) vs. 
dout for C — 2, 3, and 4. Horizontal lines mark the maxima 
provided by Eq. ©: C = 2 (red dotted), C = 3 (green dash 
dotted), and C — 4 (blue dashed). Each line intersects the 
corresponding S ■ I maxima curve at a particular d ou t, that 
coincides with the threshold of complete synchronization, as 
shown by the corresponding vertical lines in (d). 
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i.e., there exists a value of inter-community links at which 
£ reaches its maximum value (the second derivative of 
£{£out) is indeed negative). 

We now show that the maximization of £ has a dy- 
namical counterpart for the analytically treated exam- 
ple. To this purpose, we consider a modular network of 
N nodes initially arranged in two cliques c\ and C2 of 
equal size N c = f (Fig. QJa)), with ^ t = ^ 2 ut = and 
C = C = ^c(-/Vc - l)/2 such that the total number of 
links is L = N C (N C — 1). The units of the network are 
taken to be phase oscillators evolving according to the 
Kuramoto model 1141: 
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where 0j is the angular velocity of the i-th oscillator, Ui 
its natural frequency, dij is the coupling strength between 
oscillators i and j (in the following dij = di n if the nodes 
i and j are members of the same community, and d^ = 



dout otherwise) , and = 1 if i and j are connected and 
otherwise. Oscillators i = 1, . . . ,N C (i = N c + 1, . . . , JV), 
initially assigned to clique c\ (c^) and colored in blue 
(red) in Fig. [TJa-c) , have uji uniformly distributed around 
uj Ci = 0.8 (uj C2 = 0.2) within a range of ±0.1. 

To explore the extent of validity of Eq. (j3]), we start 
from such a completely segregated configuration, and 
gradually increase the ratio £ ou t/2L by i) randomly se- 
lecting a link in each module, ii) deleting it, and Hi) form- 
ing two new inter-community links by pairwise connect- 
ing the ends of the deleted links. The rewiring process 
is then repeated up to the point in which the modularity 
index Q cancels, and both modules become statistically 
indistinguishable (Fig. [ljb)). According to Eq. ([3]), the 
balance for a C = 2 module network is found at the ratio 
£ out /2L = 1 - ^Tj2. 

Furthermore, the tuning of the coupling strengths 
di n and d out allows one to control the degree of intra- 
community (cluster) and inter-community (complete) 
synchronization in each given configuration. Namely, we 



compute the order parameter R = (| l/N^jLi 
(squares in Figs. QJg-h)) and the frequency standard de- 
viation a u (circles in Figs.QJg-h)), averaged over a proper 
time window T, to account for the phase and frequency 
synchronization (in our simulations T = 200 t.u.). In 
Fig. [TJg), di n is increased up to d* n = 0.14, that con- 
stitutes the threshold for cluster synchronization in the 
case of the network structure of Fig. [IJ a). Fig.QJd) shows 
that the instantaneous frequencies 4>i of all oscillators are 
locked to their respective cluster frequencies, 0.2 (in red) 
and 0.8 (in blue). On the other hand, when the modular- 
ity index vanishes (as in the network of Fig.QJb)), we set 
din to d* n and vary d ou t to find the threshold for global 
synchronization at d* out = 1.10 (Fig. QXh)). As shown in 
Fig. QJe), for d ou t = 1.40 (sufficiently above the transi- 
tion), the network is performing a collective oscillation 
at the mean frequency of the two clusters (0.5). 

This way, functional segregated and integrated states 
are guaranteed by choosing di n and d ou t above the thresh- 
old couplings for cluster synchronization and global 
synchronization. Depending on the fraction of inter- 
community links, the competition between the dynamics 
of the individual clusters and the whole network will give 
rise to a certain degree of functional segregation and in- 
tegration. Figure [IJf) shows the instantaneous frequen- 
cies of the oscillators coupled according to the scheme of 
Fig. [TJc) for £ out /2L = 0.25. One can easily see that the 
behavior is far from being totally integrated or segregated 
but, instead, the instantaneous frequencies of the oscil- 
lators undergo high amplitude oscillations around their 
natural values. 

To quantify the degree of dynamical segregation S and 
dynamical integration /, we calculate the ensemble aver- 
ages S = (Si) and I = (Ii). Si and Ii are the dynamical 
segregation and integration measures for the individual i 
oscillators, and are defined as Si = (|<& - " C1 3 ° 2 |) T (the 
time-averaged absolute distance between the oscillator's 
instantaneous frequency and the mean frequency of the 
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two clusters), and I; = {\<t>% — £(s,, Ci)w ci —6(si, C2)w C2 \} T , 
where Si is the community of which the oscillator i is 
a member, and 5(si,Sj) the Kronecker delta function. 
The dependence of these quantities on the mixing ratio 
£ out /2L is then normalized, and reported in Fig- i) for 
din = d* n and d ou t = 1-40 > d* outl showing a monotonous 
decreasing behavior of S (squares), as the modularity 
Q (dashed line) vanishes, while a monotonous increas- 
ing trend of / (circles), up to saturation when the whole 
network is fully synchronized. 

Consequently, the dynamical segregation/integration 
trade-off, measured as the product of S and /, gives in- 
formation about the existence of a level of topological 
mixing for which functional clustering balances global 
synchronization, as shown in Fig|5Ja) (triangles). An al- 
ternative way to measure the combination of dynamical 
segregation and integration is by means of the complexity 
index E, introduced in Ref. [10| in the context of oscil- 
latory networks. Here, E is calculated using the Shan- 
non entropy of the distribution P(uj) of the average fre- 
quencies of all oscillators as E = (— YlTLi P hiP;)/ lnm, 
where m is the number of bins in the histogram of P (uj). 
By definition, E should be close to zero for narrow distri- 
butions, while it should take large values for broad dis- 
tributions reflecting the emergence of complexity. The 
index E is plotted in Fig. Ufa) (full circles) together with 
the product of the dynamical indexes S and / showing 
a noticeably similar behavior. It is important to remark 
that the calculation of S ■ I implies much simpler opera- 
tions, once the output of system Q is available. 

To more closely inspect the relation between topol- 
ogy and dynamics, we study in Fig. [2jb) the influence 
of the coupling d out on the balance between segregation 
and integration. The main observation is that the num- 
ber of inter-community links £ ou t that compromises the 
balance between the two competing synchronization pro- 
cesses decreases as the link strength d ou t increases. We 
confirm that this trend also holds for modular networks 
with C > 2. Namely, we have constructed three and 
four clique networks and performed the same analysis 
done for C = 2, obtaining the fraction £ out /2L at which 
the maxima of the S ■ I curves occur. Results are re- 
ported in Fig. [2jc) for the three values of C (triangles 
for C = 4, squares for C — 3, and circles for C = 2). 
Therefore, the functional balance depends on d out and 



arises for a given ratio £ ou t/2L where the frequency syn- 
chronization within the modules is still effective and, at 
the same time there is significant coordination between 
the modules characterized by the much richer behavior 
with the presence of several time scales (as it can be ob- 
served from Fig. HJa) where the entropy of the distribu- 
tion of frequencies exhibits a maximum). However, the 
balance between modularity and algebraic connectivity 
is maximum at one particular value of the mixing ratio 
as expressed by Eq. ^ (see Fig.^c)). Strikingly, the in- 
tersection of the analytical C line CLT/ 2 ^ = 1 - VVC 7 
with the curve of the maxima location of the S ■ I curve 
predicts the value of d ou t that perfectly coincides with 
the coupling threshold for the onset of the global syn- 
chronization, as shown by the vertical lines of Fig. EJd). 
Qualitatively similar scenarios have been observed for 3D 
chaotic oscillators, and for 2D excitable units. 

In conclusion, we introduced an easily computable 
topological measure from the knowledge of the adjacency 
matrix of the graph, and from a given partition, which 
locates the effective crossover between segregation and 
integration in a modular network, and shows that its 
maximum has a dynamical counterpart in the tradeoff be- 
tween the onsets of cluster and complete synchronization 
of networked phase oscillators. Our results can therefore 
enlighten over the performance of biological systems, that 
have to organize their structure and function to simul- 
taneously perform specialized computations at smaller 
scales and bind information at larger ones. As such, it 
can be applied, in principle, to any real world modu- 
lar network, namely to evaluate the extent to which a 
specific configuration has been optimized for information 
processing, as well as in connection with genetic, or sim- 
ulated annealing algorithms, for the generation of the 
optimal modular structure with a given number of nodes 
and links. 
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